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Abstract 

We discuss the supersymmetric 0(a s ) QCD corrections to e + e~ — ► 
Mj (i,3 = 1, 2) and to q t -> q'xf, QXk (h3 = 2 \ k = 1 • • • 4 ) within the 
Minimal Supersymmetric Standard Model. In particular we consider the 
squarks of the third generation U and bi including the left-right mixing. 
In the on-shell scheme also the mixing angle has to be renormalized. We 
use dimensional reduction (which preserves supersymmetry) and com- 
pare it with the conventional dimensional regularization. A detailed 
numerical analysis is also presented. 



1 Introduction 

In supersymmetry (SUSY) one has two types of scalar quarks (squarks), q L 
and q R , corresponding to the left and right hclicity states of a quark. q L and 
q R , however, mix due to the Yukawa coupling to the Higgs bosons, which is 
proportional to the mass of the quark. One therefore expects large mixing in 
the case of the stop quarks so that one mass eigenstate (m^ ) might be rather 
light and even reachable at present colliders. The sbottoms &£,, b R may also 
strongly mix for large tan/3. The mass matrix in the basis (q L , q R ) is given 
by: 

M * = ( m i +m t +D * ^-"{Zl})) (1) 

V^O** m\ b + ml + D 2 ) ' 

Here ttLq, m (j, m £)i an d A q are SUSY soft -breaking parameters, fi is the 
Higgsino mass parameter, and tan/3 = j 2 -- D 1 and D 2 are the D terms: 
D 1 = TO^cos2/3(/g L - e q sin 2 8 W ), D 2 = m| cos 2/3 e q sin 2 d w , with i| L the 
third component of the weak isospin of q. In the off-diagonal elements of Eq. ([!]) 
cot j3 enters in the case of the stops and tan/3 in that of the sbottoms. Diag- 
onalizing the matrix one gets the mass eigenstates = q L cos 9q + q R sin 9q, 
q 2 = —q L sinflq + q R cos9q with the masses m^, rriq 2 (with < rriq 2 ) and 
the mixing angle Oq. 

f Talk presented at the International Workshop on Quantum Effects in the MSSM, Septem- 
ber 9 - 13, 1997, Barcelona, Spain. 
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Conventional QCD corrections to squark pair production e + e~ — > q^j 
= 1, 2) can be very largely. The SUSY-QCD corrections including squark 
and gluino exchange will be discussed here following closely ref . □. These cor- 
rections were also treated in 13. The new feature in the calculation of the 
SUSY-QCD corrections is that in the on-shell scheme a suitable renormaliza- 
tion condition has to be found for the mixing angle 9 q because the tree-level 
amplitude explicitly depends on it. We will explain this in detail below. 

The SUSY-QCD corrections to the squark decays into chargino or m 
tralino, g s — > q'xf, qXk (*>J = 1) 2; k = 1 . . .4), have been calculated in 
and will also be discussed in the following. Here the dependence on the na- 
ture of the charginos/neutralinos (gaugino-like or higgsino-like) is particularly 
interesting. 

We work in the on-shell scheme and use dimensional reduction (DR) to 
regularize the integrals, which is necessary to preserve supersymmetry (at least 
up to two loops). We will comment on the differences between this and the 
dimensional regularization scheme used in the Standard Model. 

2 The Production Process e + e — + q^j 

The cross section at tree level is given by: 

cr° (e + e~ -> ftfj.) = —A-/ 2 [e 2 q <Sy - T lZ <,,",,*,, + T zz o|] (2) 

With rp _ Ve s{s-m%) 

lZ ~ 8c 2 w s 2 w [{s-m%Y+Tlm%y 

zz " 256 s^ w (s - mlf + T|m| ' W 

Here Ay = (1 — /if — fi 2 ) 2 — 4/xf /Lt^ with fj^j = m|. Js. e q is the charge of the 
squarks (in units of e), v e = — 1 + 4s^, a e = —1, % = sin9w, cw — cos Ow- 
dij are the relevant parts of the couplings Zq^q*: 

a n = 4 (7 3L cos 2 9 q - s 2 w e q ) , a 22 = 4 (l| L sin 2 9 q - s 2 w e q ) , 

aii = 021 = -2/ g 3L sin2^. (5) 

The SUSY-QCD corrections in 0(a s ) consist of the conventional QCD 
corrections □ due to gluon exchange and real gluon radiaton, as well as of the 
corrections due to the exchange of a gluino and squarks, see Fig. 1. Our input 
parameters are the physical masses m qi , m q2 , rriq, m q , and the mixing angle 
9 q . We use the on-shell sheme where the masses are fixed by the respective 
poles of the propagators. 
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Figure 1: Feynman diagrams for the lowest order SUSY-QCD corrections to e+e~ — > q^j 
with squarks and gluinos in the loop. Note that there are also the corresponding diagrams 
to b), and d) for the antisquark <jj. 



In renormalizing the lagrangian we follow the usual procedure: 

£ = C + SC (6) 

with 

C = -ee q 8 l3 A» q*(id^ q } - y q 3 . (7) 

Cq, the bare lagrangian, has the same form with the bare quantities: 

e q$ij = e q<>ij + (Se q )ij , (8) 

a°- = ciij + Saij , (9) 

q*° = {l + \5Z tl )q* +5Z ll ,q* l , i ^ i' , (10) 

g° = (l + ^Z n )q l+ 5Z^q^ jjtf. (11) 

Notice that because of 0? = 6q + 59 qi Sa^ is a function of S9 q . The total 
correction in 0(a s ) can be written as: 

Aa y = ^j'+^+Ja'?, (12) 
(Ae,)« - (Se q )\f + (Se g )\f, (13) 
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where (v) denotes the vertex corrections and (w) the wave-function corrections. 

The contributions come from gluon, gluino, and squark exchange. 5a\j is due 
to the shift from the bare to .the on-shell couplings. As already mentioned, 
we use dimensional reduction □ instead of dimensional regularization. Up to 
first order this is achieved technically by taking D = (4 — r e) with r — > (see 
section 4). In the case of e + e~ — > there is, however, no difference between 
the two schemes as will be explained later. 

Let us first discuss the vertex corrections Se\j and tffljj coming from 
the exchange of SUSY particles. The gluino contribution due to the graph in 
Fig. la is given by: 

*°«* > = ~{^m~ g m q v q Sl j {2C± + C%) 

+ v q 8 l3 [(2ml + 2m 2 q + m? + m? ) C± + 2m? C% + B°(s, m 2 q ,m 2 q )] 
+ a q A% [(2m? - 2m 2 , + m? + m| ) C± + (m? - m? ) 

+2m 2 C° J +B ( S ,m 2 ,m 2 q )]} (14) 

and 

6(e q )%'*> = ^ e ,{2mgm g 4(2C+ + Cf i ) (15) 
+ * y [(2m? + 2m 2 , + m | + m | ) C + + 2mgC° + m 2 ,, m 2 ,)]} 

with v q = 2I* L - As 2 w e q , a q = 21 , Sf x = - sm29 q = -Sf 2 = A\ 2 = A\ x , and 
S\ 2 = S 21 — — cos2#,j = — A\ x = A 22 . The functions are defined by 

C+ = ^, <T-<^. ,16, 

B° and Crj 1 ' 2 are the usual two- and three-point functions as given, for in- 
stance, intl The arguments of all C-functions are (to?, s, to?. , to?, to 2 ,, to 2 .). 
The squark exchange graph Fig. lc is proportional to the four-momentum of 
Z° , and therefore does not contribute to the physical matrix element. The 
wave-function corrections (Figs. 1 b ,d) can be written as, using Eqs. (^) to 
(0) (i^i'Jjtf): 

Sa^ = \{5Za + 8Zjj)a,ij + 5Z l > i a V] + SZjij ay . (17) 

An analogous formula holds for S(e q )^ with — > e q S^j. 
One obtains from Fig. lb: 

= -Refl^KD + Ef^m?.)]^, 



ggK) a , g£K) a i , ltn 

Tit — mi mi — mi J 



and 



S(e 9 t' S) = -e 9 Re{s:f(m|)}, (19) 
- o 6 _ 9 o Re | sffi (m| a ) - sffi (mf x ) } , (20) 

where E^(m 2 ) are self-energies and T,'^\m 2 ) — dY±\f'(p 2 ) / 'dp 2 \ p i =m 2 . Notice 

that S(e q )\j = because the contributions with the squark loop attached 
at either external squark line in Fig. Id cancel each other. The wave-function 
correction Sa[j'^ due to Fig. Id plays an important role in the renormalization 
of the squark mixing angle 6q. 

2.1 Renormalization of the Mixing Angle 9 q 

The total correction Aaij , Eq. ( |l2| ) , using Eq. ( |l7|) can be written as (i ^ i', j ^ 
f) 

Aa^ = Sa\y + h(5Z ti + SZjj) + 5Z Vi a Vj + 5Zj>j a ijt + 5af) . (21) 
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Notice that the first part of the right-hand side, Sa^ + \(bZa + 8Z 
is already free of ultra-violet divergencies. Hence, the second part of Eq. (21) 
has to be finite, too. We therefore may require for i = 1 and j = 2 

Sa\ 2 ' = (a 22 - an)S9q = - (5Z 21 a 22 + SZ 12 an) . (22) 

One can easily see that Aaij is then also finite for all i, j. The condition, 
Eq. ([22]), means that the non-diagonal self-energy graphs Fig. lb and Id can- 

(f)\ j — V 

eel the counterterm Sa 12 in Eq. (21). Notice also that the total squark con- 



tribution Aafj is zero. Other authors used the same basic idea but took, 

for instance, the condition analogous to Eq. d22|) for Sa^ or 8a 2 ® 2 \ see@, or a 
similar condition valid at a point Q 2 , see ref. tl The differences between these 
schemes are, however, numerically very small. 

2.2 Total QCD Correction in 0(a s ) 

The total QCD correction Act to the cross section is 

Aa = A<t (s > + Aa^ , (23) 
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as Act® = in our renormalization scheme of the squark mixing angle. The 
gluon contribution factorizes: 



4 as 

3 TT ' 



(24) 



where Ay is given in ref. □. The total gluino contribution is given by: 



s J 



|2e,(Ae g ) 



ij 



2T ZZ aijAa tj 



(g) 



-T jZ [e q S tJ Aaf + (Ae q )fa tJ ]} (25) 



with 



(§) 



-Re{i[S^(m|)+E^(m|)]a 



4 a s m,gm q 



B\m\,ml, m 2 q ) [(-1) 4+1 1a iv cos 2% - a vv sin20 ? ~] 
+ B^m^TO^m^aii sin 20 ? } (26) 



(i^t') and A(e,)<?=(a e ,)g> + (*e,)#>. 



2.3 Discussion 

First, we have calculated the SUSY-QCD corrections to the cross section of 
e+e~ — > iiti in the LEP energy range \/s < 200 GcV. We have found that, 
whereas the conventional QCD correction may be rather large, the gluino cor- 
rection is only about 1% of the tree-level cross section, quite independent of 
m$ . The correction due to gluino exchange is, however, not negligible (2 — 8%) 
in the energy range of a linear e + e~ collider (y/s = 500 — 2000 GeV). 

The -^"dependence f th e SUSY-QCD corrections to the cross section 
a(e + e- -> tik) is shown in Fig. 2 for m fl = 100 GeV, m ?2 = 400 GeV, 
■nig = 300 GeV, and cosOi = l/\/2. The peak at \/s = 350 GeV is due to the 
ti threshold. In Fig. 3 we show the cos 0^ dependence of the corrections for this 
process at y/s = 500 GeV for the same masses of the stops and the gluino as 
in Fig. 2. Whereas the gluon correction has the same behaviour in Of as the 
tree-level cross section, the gluino correction is different. In Figs. 4 and 5 we 
exhibit the corrections to cr(e + e~ — > t^) and a(e + e~ — > i2i2), respectively, at 
y/s = 2 TeV for m h = 400 GeV, m t ~ 2 = 800 GeV, m g = 600 GeV. The gluino 
contributions can go up to about —10%. Fig. 6 shows the dependence on the 
gluino mass. It is interesting to notice that the gluino correction decreases 
very slowly with the gluino mass. 
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Figure 2: SUSY-QCD corrections 

Sa g/ a tree and fog^tree for e + e - _> ^ 

as a function of for cos 9^ = 1/ \/2, 
m fi = 100 GeV, m t ~ 2 = 400 GeV, and 
m s = 300 GeV. 
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Figure 3: SUSY-QCD corrections 

Sa g/ a tree an( j fog^tree for e+e - _> ^ 

as a function of cos 0j for = 500 GeV, 
mfi = 100 GeV, m t ~ 2 = 400 GeV, and 
m- a = 300 GeV. 
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Figure 4: SUSY-QCD corrections Sa 3 and 
Scr 3 as a function of cos 0j for e+e - — » £1*2, 
= 2 TeV, m fl = 400 GeV, m ( ~ 2 = 
800 GeV, and m 5 = 600 GeV. 




Figure 5: SUSY-QCD corrections 

8a g/ a tree and fog j a tree for £ + £ - _> ^ 

as a function of cos 0{ for v 7 ! = 2 TeV, 
m fl = 400 GeV, m ?2 = 800 GeV, and 
m s = 600 GeV. 
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Figure 6: Dependence of the SUSY-QCD 
corrections 8<r 9 /a tree and 8a 9+g /a tree on 
the gluino mass for e+e - — > t\t\, for 
= 500 GeV, m fl = 100 GeV, = 
400 GeV, cos 6 { = l/s/2. 



3 Squark Decays into Charginos and Neutralinos 

In the following we discuss the SUSY-QCD corrections for the decays: 

U -» bxf, k -» 4*7, (27) 
*< -> *Xfc, k -» (28) 

with i, j = 1,2 and fc = 1 . . .4. The supersymmetxic QCD corrections were 
calculated for m g = and Xi being a photino in ref.E2l, and taking into account 
squark mixing, quark masses (i.e. Yukawa couplings), and general gaugino- 
higgsino mixing of charginos and neutralinos in refs.u andu. The decay width 
at tree-level for tj — > bx^ is given by: 

2/2 2 2 \ 

T°(ti^bxi) = ' — 



167tto? 



• ([(4) 2 + (4) 2 ] X - 1 '<%'''<<'' v) (29) 
- ijjj — m? + and k(x, y, z) = \(x — y — z) 2 — Ayz] 1 / 2 . The t*-b~xt 



couplings £jf- and k\j read, for instance, for t\ — » bjj7 • 

4 = -^iCos%+ mt . K-2sinfl f , (30) 
V^TOvvsinp 



fefj = -jJ^ -U^cosd;, (31) 

V 2 tow cos p 



where U and— are the matrices diagonalizing the charged gaugino-higgsino 
mass matrixO. 
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The 0(a s ) SUSY-QCD corrected decay width can be written as: 

r = r° + <5r (l,) + 5r (w) + sr {c) + <5r (real gluon ) , (32) 

where the superscript v again denotes the vertex correction (Figs. 7a,b) and w 
the wave-function correction (Figs. 7c-g). ST 1 -^ corresponds to the shift from 
the bare to the on-shell couplings, taking into account the renormalization of 
the quark mass and the squark mixing angle. (5r( real s luon ) [ s the correction due 
to real gluon bremsstrahlung and cancels the infrared divergencies. 
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Figure 7: Vertex and wave-function corrections to squark decays into charginos and neutrali- 
nos. 



The procedure of the calculation is completely analogous to that discussed just 
before in section 2. The complete formulae for the different correction parts in 
Eq. (|3^) are given in ref. We want to note that, contrary to the production 
process e + e~ — > q^j, the corrections to the decay widths of % — ► q'xf an d 
Qi ~* IXk arc different in the dimensional regularization and in the dimensional 
reduction scheme. (At first order the difference is finite.) This is because of 
the quark wave-function correction due to gluon exchange, Fig. 7c. The quark 
self-energy corresponding to Fig. 7c is given by: 

n (s) (fc 2 ) = 2L [ 21/zB 1 + 2(jS - 2m q )B° - r(ft - 2m,)] (33) 
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with B n = B n (k 2 , A 2 , to 2 ,) and the gluon mass A — > 0. This leads to the quark 
wave-function rcnormalization constants due to gluon exchange 



3 7T 



B° + B 1 -2m 2 q (B° -B 1 )- - 



(34) 



with B n = £™(to 2 ,A 2 ,to 2 ), B n = B n (to 2 , A 2 , to 2 ). SZ l and SZ R arc defined 
by the usual relation between the unrenormalized quark field q° and the renor- 
malized one, q° = (1 + \&Z L P L + ^SZ R P R ) q. Note the dependence on r in 
Eqs. ( j33|) and ([mJ), where r = in the dimensional reduction and r = 1 in 
the dimensional regularization scheme. Note, however, that there is no such 
difference for the squark self-energy graph due to gluon exchange. 



3.1 Numerical Results 

Let us first discuss the decay t\ — > bxt ■> where we take to~+ = 100 GeV, 
tan/3 = 2, = 600 GeV, = 450 GeV, = 470 GeV, and cos% = 
-0.9. We study three cases: M < |/i| (M = 95 GeV, fi = -800 GeV), 
M ~ (M = 100 GeV, fi = -100 GeV), and M > \n\ (M = 300 GeV, 
/i = -89 GeV). We use the GUT relations: M' ~ 0.5 M, m~ g ~ 3.5 M. 

In Fig. 8 the dependence of the SUSY-QCD corrections on the stop mass is 
exhibited for cos 9^ = 0.6. Notice the pronounced dependence on the nature of 
the chargino. The corrections are largest (~ —25%), if the chargino is higgsino- 
like (|/i| -C M) due to the large top Yukawa coupling. If xt is gaugino-like 
(M <C \n\) the corrections are between +20% and —10%. 

In Fig. 9 we show the SUSY-QCD corrected widths together with the tree- 
level widths as a function of cos 9^ for to^ = 200 GeV and the other parameters 
as in Fig. 8. Again, the corrections are biggest in the case of a higgsino-like 
chargino. The behaviour of the cos6>j dependence reflects the fact that if 
ti ~ in (cos#j ~ 0) it strongly couples to the higgsino component of xt , and 
if t\ ~ tL {cos9i ~ ±1) it strongly couples to the gaugino component. 

In Fig. 10 we show <5r/r° [%] as a function of for t\ — * tx±, taking 
to^o = 80 GeV, tan/3 = 2, — 600 GeV, and cos% = 0.6. Again we observe 
that if Xi is higgsino-like (|/z| <C M) the corrections are about —20%. 

We have also studied the dependence on the gluino mass. In Fig. 11 wc 
show a plot where <5r/r° is exhibited for t\ — > bxt an d t% — > tx% as a function of 
mg for m fl = 300 GeV, cos% = 0.6, tan/3 = 2, and ^ = -100 and -800 GeV. 
M is fixed by M ~ 0.3 mg. Notice that the SUSY-QCD corrections are still 
important for nig ~ 1 TeV and no decoupling of the gluino mass can be seen. 
This is also the case if we relax the condition M ~ 0.3 TOg and keep the chargino 
(neutralino) mass fixed. 
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Figure 8: SUSY-QCD corrections to the 
width of tj — + b\f as a function of , for 
m.+ = 100 GeV, cos6» t - = 0.6, tan/3 = 2, 
and various (M, fi) [GeV] values. 
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Figure 9: Tree-level (dashed lines) and 
SUSY-QCD corrected (solid lines) decay 
widths of — » as a function of cos 9^, 



for mr 



200 GeV, m.+ 



100 GeV, 



tan/3 = 2, and various (M, /i) [GeV] values. 




<5r/r 



-20 



(300, -85) 



300 



400 500 

m fi [GeV] 

Figure 10: SUSY-QCD corrections to the 
width of t 1 — + tXi as a function of mj , for 
m-o = 80 GeV, cos0 t ~ = 0.6, tan/3 = 2, and 
various (M, £t) [GeV] values. 
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Figure 11: SUSY-QCD corrections to the 
widths of tj — > 6x)^ (solid lines) and t x — » 
txj (dash-dotted lines) as a function of m,g, 
for m fl = 300 GeV, cos% = 0.6, tan/3 = 2, 
M ~ 0.3m s . 
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4 Dimensional Reduction Technique 

The regularization by dimensional reduction was proposed byi. It means that 
only the space-time dimensions (the coordinates a; M and momenta p M ) are con- 
tinued to D = 4 — e dimensions, whereas the vector fields and spinors remain 
four-dimensional. Following it is convenient to write the four-dimensional 
vector field as V M = (Vj, V a ), where Vi is a D-dimensional vector, and 
V a is e-dimensional behaving as e scalars. Moreover, one has 7 M = (7 1 ,7 <T ). 
Note that a; M = (x l , 0), <9 M = (<9 l , 0), and = (p l , 0). As a consequence the la- 
grangian C can be decomposed as C = C^+C^ e \ where D- D > is the lagrangian 
of the conventional dimensional regularization. Therefore, to each interaction 
term of a vector field there is a corresponding "e scalar" interaction term, ex- 
cept for the vector-scalar-scalar interaction because V 11 cf>* <9 M 4> = V 1 <fi* di <fi 
(and no e term), with <p being a scalar field. Therefore, in this case there is 
no difference between dimensional regularization and dimensional reduction. 
There is, however, a difference in the case of the interaction of a fermion with 
a vector field. For instance, the fermion self-energy, Fig. 7c, receives a contri- 
bution due to e scalars in the loop of §j(^ — 2rn 9 ). This is just the expression 
which cancels the i — dependent term in Eq. (|3^) for r = 1 in order to get the 
result of dimensional reduction (r = 0). Thus at the one-loop level the "e- 
scalar" technique is equivalent to performing the algebra in the numerator of 
the integrand in four dimensions and making the integration in D dimensions, 
or equivalently taking D = 4 — re with r — + 0, as we did in our calculations. 
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